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ABSTRACT 


The  finite-difference  treatment  of  the  elastoplastic 
problem  of  a  thick-walled  cylinder  is  nor., -all  -  based  on  the 
differential  equations  for  the  displacemen'  actor;  hence, 
determination  of  stresses  and  strains  requires  numerical 
differentiation.  However,  good  results  in  differentiation 
are  not  provided  by  the  computer  unless  a  rather  fine  grid 
is  used.  With  the  alternative  method,  developed  in  this 
report,  incremental  stresses  and  strains  at  each  nodal  point 
are  directly  used  as  variables;  hence,  numerical  differen¬ 
tiation  in  the  evaluation  of  stresses  and  strains  is  not 
required.  An  effective  incremental  theory  has  been  developed 
for  solving  the  problem  of  partially  yielded,  thick-walled 
cylinders  made  of  strai n-ha* deni ng  materials  subjected  to 
any  combination  of  internal  pressure,  external  pressure,  and 
end  loading.  The  theory  developed  in  the  Research  Directorate 
U.  S.  Army  Weapons  Command,  includes  Prandtl -Reuss '  incre¬ 
mental  stress-strain  laws,  Von  Mises'  flow  criterion,  and 
the  strain-hardening  and  compressibility  properties  of  a 
material.  Since  the  consideration  of  stress  and  strain 
history  is  involved  in  the  analysis,  the  present  theory  is 
particularly  suitable  for  predicting  stress  and  strain  dis¬ 
tribution,  and  location  of  the  elastic-plastic  boundary  of 
a  thick-walled  cylinder  subjected  to  nonproportionate  loading. 
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NOMENCLATURE 


a  ,b 

A 

E 


Inside  and  outside  radius  ot  a  thick-walled 
cylinder,  respectively. 

Point  on  curves  indicating  initiation  of 
yielding  in  the  thick-walled  cylinders. 

Young's  modulus. 


F 


P*.P. 


Difference  between  axial  load  transmitted 
through  cross  section  of  a  cylinder  and  P ^ rra 2 

Internal  and  external  pressure,  respectively. 


r,p  =  Radius  and  dimensionless  radius,  respectively 

a  =  Normal  stress. 


a..  =  Stress  tensor. 

■  J 

•' ;’°2,a3  =  Principal  stresses. 


a 


c 


a 


G 


e 


Ci  j 


e  ,e  ,e 

l  i  3 


r  •  -  " 
f.  ,  6. 


Effective  stress. 

Yielding  stress  in  tension. 

Yielding  strain. 

Normal  strain. 

Strain  tensor. 

Principal  strains. 

Elastic  and  inelastic  normal  strain, 
respectively. 


Effective  strain. 


Strain-hardening  factor  for  material. 


Poisson's  ratio. 

boundary^  l'SS  rad’US  °f  ^  1  a  s  t  i  c  -  pi  a  s  t  i  c 

respectively. '^reSS  "ormal  s‘ra1n, 

The  deviator  components  of  strain  in  the 
direction  of  r,  e,  and  z,  respectively 

I?®j!rat0r/0mponents  of  stress  in  the 
irections  of  r ,  8,  and  z,  respecti vely . 


INTRODUCTION 

Preliminary  Statement 

When  a  thick-walled  cylinder  is  deformed  under  c ont  >  n- 
uously  increasing  pressure  from  an  initial  value  ze~o,  the 
entire  thick-walled  cylinder  is  at  first  elastic,  and  the 
stresses  and  strains  at  any  point  in  the  hollow  cylinder  a^e 
given  by  Lame's  equations.1^’3  These  equations  are  val’d 
until  the  pressure  reaches  the  critical  value  at  which  the 
material  at  the  bore  of  the  thick-walled  cylinder  becomes 
inelastic.  As  the  pressure  increases  above  this  critical 
value,  the  inelastic  region  spreads  outward  from  the  bore 
through  the  thick-walled  cylinder  A  section  of  a  parfaMy 
yielded  hollow  cylinder  in  which  the  inelastic  region  has 
spread  to  a  radius  c  is  shown  in  Figure  1  As  the  pressure 
is  increased,  the  elastic-plastic  boundary  is  continuously 
moved  outward  until  it  coincides  with  the  outer  surface;  the 
hollow  cylinder  is  then  said  to  be  comp’etely  yielded. 

The  fact  that  the  internal  pressure  in  a  thick-walked 
cylinder  can  rise  considerably  above  the  va’ue  for  initial 
yielding  before  the  thick-walled  cylinder  becomes  completely 
inelastic  is  of  practical  interest  for  design  of  high  pres¬ 
sure  tubes  and  gun  barrels  If  a  thick-walked  cy1inder  is 
expanded  so  that  the  region  of  yielding  penetrates  to  a  cer¬ 
tain  depth  and  the  internal  pressure  is  then  removed,  desir- 
able  residual  stresses  may  be  produced.  For  certain  residual 
stress  distributions,  if  the  inter, , a  1  pressure  ’S  reapplied, 
the  thick-walled  cylinder  will  not  yield  again  until  a  pres¬ 
sure  equal  to  the  maximum  expanding  pressure  previously  ap¬ 
plied  has  bee.  reached.  The  elastic  strength  of  a  partially 
yielded  thick-walled  cylinder  thus  may  be  considerably  greater 
than  that  of  an  untreated  thick-walled  cylinder  of  the  same 
dimensions.  Hence,  autofrettage  of  thick-wa1 ’ed  cylinders 
is  becoming  a  popular  industrial  orcctice,  pa  i  cut  ar  )y  for 
gun  barrels. 

In  almost  all  theories  presented  in  the  literature,  the 
thick-wailed  cylinders  are  assumed  to  be  made  of  nonstrj’r>. 
hardening  materials  Hill,  tee,  and  Tupper  presented  com¬ 
pressible  theories  for  tnick  tubes  under  plain-strain  assump¬ 
tion"  and  for  closed-end  tb' ck-wal led  cylinders*  made  of 
nonstrain-hardening  materials.  The  Tresca  yield  condition 
and  the  Prandt 1 -Ruess  stress-strain  relations  were  useJ  in 
the'r  investigations. 
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FIGURE  1  Partially  Yielded  Thick-Walled  Cylinder 
Subjected  to  Internal  &  External  Pressures 
and  End  Loads 


Hodge  and  White6  proposed  compressible  theories  for 
thick-walled  cylinders  made  of  nonstrain-hardening  materials 
with  zero  axial  strain.  Von  Mises'  yield  criterion  was  used 
in  their  theories.  They  found  little  difference  between  the 
incremental  theory  based  on  the  Prandtl -Reuss  stress-strain 
relations  and  the  theory  based  on  Hencky's  total  stress-st^a-* n 
relations . 

MacGregor,  Coffin,  and  Fisher7  have  presented  a  com¬ 
pressible  theory  for  open-end,  thick-walled  cylinders  maoe 
of  nonstrain-hardening  materials.  Von  Mises'  yield  criterion 
and  Hencky's  total  stress-stra in  relations  were  assumed  to 
be  valid  in  their  investigation. 

Or  the  basis  of  the  assumption  that  the  axi.al  strain  is 
purely  elastic,  that  the  axial  stress  can  be  expressed  in 
terms  of  the  axial  strain,  and  that  the  radial  and  circumfer¬ 
ential  stresses  can  be  determined  by  means  of  Hooke's  law, 
Bland8  presented  a  theory  for  the  thick-walled  cylinders 
made  of  linear  strain-hardening  materials.  He  used  Tresca's 
yield  criterion.  In  his  investigation,  the  radial  stress  is 
found  to  be  independent  of  axial  strain.  Also  with  this 
theory,  the  pressure  necessary  to  locate  the  elastic-plastic 
boundary  at  a  given  radius  can  be  predicted.  (The  elastic- 
plastic  boundary  is  independent  of  end  conditions.) 

Steele9  proposed  a  closed-solution  theory  for  open-  and 
closed-end  thick-walled  cylinders  made  of  linear  strain-hard¬ 
ening  materials.  He  used  Tresca's  yield  criterion  and  the 
Hencky's  total  stress-strain  relations.  According  to  his 
theory,  discontinuities  have  been  found  ut  the  elastic-plastic 
boundary  for  axial  stress  and  bore  circumferential  strain  if 
the  elastic  material  and  plastic  material  are  assumed  to  be 
compressible  and  incompressible,  respectively 

Smith  and  Sidebottom10  have  presented  an  incompressible 
theory  for  thick-walled  cylinders  made  of  linear  strain-hard¬ 
ening  materials  and  subjected  to  internal  pressure,  external 
pressure,  and  axial  load.  They  used  Von  Mises'  yield  con- 
ditior  and  Hencky's  total  stress-strain  relations.  Their 
theory  can  be  corrected  to  predict  pressures,  stresses,  and 
strains  for  thick-walled  cylinders  made  of  compressible  ma¬ 
terials  as  long  as  the  elastic-plastic  boundary  has  not 
reached  the  outer  radius.  Experimental  data11  have  shown 
that  the  external  strains  have  been  accurately  predicted  by 
this  theory. 
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To  the  author's  knowledge,  neither  a  closed-form  com¬ 
pressible  solution  nor  a  numerical  compressible  solution  has 
been  proposed  to  solve  the  problem  of  stress  and  strains  in 
thick-walled  cylinders.  These  cylinders  are  made  of  strain- 
hardening  material  and  subjected  to  nonproporti onate  loading 
for  either  open-end  or  closed-end  conditions  based  on  Prandtl- 
Reuss  incremental  stress-strain  relations  and  Von  Mises'  yield 
condition.  The  aim  of  this  investigation  was  to  derive  a  gen¬ 
eral  theory  for  thick-walled  cylinders  satisfying  all  the  re¬ 
quirements  listed  above. 

Stress-Strain  Relations 

The  most  distinguishing  characteristic  between  elastic 
and  inelastic  analyses  is  the  loading  path  dependence.  If  a 
member  is  loaded  elastically,  the  final  state  of  stress  and 
strain  is  independent  of  the  loading  history,  and  the  state 
of  strain  is  always  the  same  for  the  same  state  of  stress  and 
vice  versa.  Stated  another  way,  in  the  elastic  range,  strains 
or  stresses  are  uniquely  determined  by  the  stresses  or  the 
strains.  However,  for  materials  in  the  inelastic  range,  the 
final  state  of  stress  cannot  be  determined  by  the  final  state 
of  strain,  or  vice  versa.  Therefore,  in  the  inelastic  range, 
strains  are  in  general  not  uniquely  determined  by  the  stresses, 
or  vice  versa,  but  depend  on  the  history  of  loading  or  how  the 
stress  state  was  reached. 

The  inelastic  stress-strain  relations  for  a  given  material 
are  empirically  determined.  Usually  tension  and  compression 
specimens  are  tested  to  determine  the  stress-strain  relation 
for  a  uniaxial  state  of  stress. 

In  Figure  2,  curve  OAB  represents  the  constitutive  re¬ 
lation  of  a  member  made  of  strain-hardening  material  and  sub¬ 
jected  to  uniaxial  state  of  stress  under  the  condition  that 
the  loads  on  the  member  cause  the  stress  in  all  inelastical  ly 
deformed  volume  elements  to  be  monotonically  increasing  in 
magnitude.  From  0  to  A,  the  behavior  of  the  material  is  ap¬ 
proximately  linearly  elastic.  The  stress  at  A  may  be  con¬ 
sidered  the  yield  stress,  ae,  since  from  A  to  any  point  B  on 
the  curve  the  behavior  of  the  material  is  inelastic.  In 
general,  the  loading  function  for  a  given  material  for  other 
states  of  stress  must  be  predicted  from  the  loading  function 
for  uniaxial  state  of  stress. 
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FIGURE  2  Stress-Strain  Diagram 


For  multiaxial  state  of  stress,  at  each  point  in  a  given 
member  and  at  each  stage  of  plastic  deformation,  a  numDer, 
k>0,  and  a  function, f(o. .) ,  exist  such  that  further  inelastic 
deformation  takes  place1,3  only  if 


f(a1j)  >k 


(1) 


A  yield  surface  in  stress  space  is  determined  by  the  condition 


f (°1 j)  -  k  (2) 

On  or  within  this  stress  surface,  the  behavior  of  the  material 
is  linear  elastic  only.  Equation  2  with  specified  function  f 
and  specified  value  k  represents  one  of  an  infinite  number  of 
yield  surfaces  for  a  given  strain-hardening  material,  depend¬ 
ing  upon  the  amount  of  Inelastic  deformation  as  well  as  on 
prior  plastic  strain  history.  The  only  one  of  these  Infinite 
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number  of  possible  relations  for  Equation  2,  generally  assumed 
to  be  known  for  a  given  material,  is  the  one  that  represents 
the  yield  surface  of  the  virgin  material.  In  general,  the 
yield  surfaces  and  resulting  stress-strain  relations  will  de¬ 
pend  upon  the  loading  path.  Obviously,  an  equation  of  state 
does  not  exist  in  the  theory  of  inelasticity.  If  the  equa¬ 
tion  of  state  did  exist,  it  would  depend  upon  infinitely 
many  state  variables.  Hence,  in  the  strict  sense,  prediction 
of  the  plastic  strain  components  e j  .  that  will  result  from  a 
given  stress  history  is  impossible. J  The  experiment  itself 
must  be  run  to  determine  e^. 

Linearly  elastic  theories  require  that  material  for  a 
given  state  of  stress  does  not  exceed  the  yield  for  that  state 
of  stress.  Inelastic  theories  require  that  the  yield  stress 
for  a  given  state  of  stress  be  known  to  determine  the  elastic- 
plastic  boundary.  Thus,  a  yield  condition  must  be  specified 
to  predict  the  yield  stress  for  a  given  state  of  stress  from 
the  known  yield  stress  obtained  from  tension  and  compression 
specimens  or  from  hollow  torsion  specimens. 

Experiments12*13  on  metals  have  shown  that  initial  yield¬ 
ing  and  subsequent  plastic  flow  are  not  affected  by  a  moderate 
hydrostatic  pressure.  This  reduces  the  dependence  of  the  yield 
condition  of  an  Isotropic  material  on  its  second  and  third 
stress  deviator  Invariants  only.16  The  two  most  accepted  yield 
conditions  are  Tresca's  maximum-shearing-stress  yield  condition 
and  Von  Mises*  yield  condition. 

Von  Mises*  yield  condition  states  that  the  second  stress 
deviator  invariant  is  a  constant,  i.e., 


—  [(o  -o  )2  +  (o  -o  )*  +  (o  -o  )i  +  6  (t  2+t  2  +  r  2)]5* 

V  y  Wj,/  V  2  yJX,  W  V^y  ly2  ^  /j 


The  constant  as  determined  from  a  tension  test  Is  o  //J. 

f? 
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Tresca's  yield  condition  states  that  yielding  occurs 
at  a  point  in  a  stressed  volume-element  whenever  the  max¬ 
imum  shearing  stress  at  that  point  reaches  a  certain  crit 
i cal  val ue ,  i . e.  , 


°max  ‘  amin  *  k2 


in  which  k  is  a  constant  usually  determined  from  a  simple 
test.  Since  Equation  4  should  be  the  same  for  all  states 
of  stresses,  k2  *  afi  where  aQ  is  the  yield  stress  in  tension. 

In  general,  the  combination  of  the  principal  stress 
that  gives  the  maximum  shearing  stress  is  not  known  before¬ 
hand.  The  condition  expressed  in  all  generality  is  given  by 


C(a1-a2)2-k22][(a2-oj)2-k22][(o3-ai)'-k22]  =  0  (5) 

Considerable  variance  Is  present  concerning  the  validity 
of  any  one  yield  condition;  however,  experimental  results11* 
Indicate  that  test  data  for  ductile  metals  have  shown  reason¬ 
able  agreement  with  both  the  Tresca's  and  the  Von  Mlses'  yield 
conditions,  these  data  are  more  favorable  to  Von  Mlses' 
yield  condition.  In  this  Investigation,  Von  Mlses*  yield 
condition  is  assumed  to  be  valid. 

The  most  popular  loading  function  for  strain-hardening 
materials  may  be  written  as  for  all  states  of  stress, 

where  the  effective  stress  <T  and  effective  strain  e  are  de¬ 
fined  by 

°*^2‘S/(ox-<’y)*  +  (Oy-o2)*  +  (o£-ox)*+6(Txy2+TyZ  *+tzx2)  (6) 


and 

(ex-€y),  +  (ey-ez Yy2 * +VZX 1 )  ( 7 ) 
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Another  loading  function  proposed  by  Ludwik18  states 
that  the  maximum  shearing  stress  as  a  function  of  the  max¬ 
imum  shearing  strain  describes  the  subsequent  flow  of  a 
material  in  the  strain-hardening  region. 

Several  theories  have  been  proposed  to  describe  the 
stress-strain  relations  of  a  volume  element  in  the  inelas¬ 
tic  range.  Two  of  these  theories,  Hencky's  total  strain 
theory  and  Prandtl -Reuss '  incremental  strain  theory  have 
received  wide  attention. 

Hencky's  stress-strain  relations  provide  a  logical  ex¬ 
tension  of  Hooke's  elastic  equations  and  embody  three  hy¬ 
potheses  : 10 

1.  Principal  axes  of  stress  and  strain  coincide. 

2.  Mohr's  circle  diagrams  of  stress  and  strain  are 
similar  at  any  stage  of  the  inelastic  deformation. 

3.  Volume  changes  are  elastic. 

In  cylindrical  coordinates,  the  first  two  hypotheses 
are  satisfied  by  the  relation 


e  -e 
r 


o  -o„ 
r  e 


e  _  ee“ez  .  ez"er  .  Yre  _  Yez  _ 

-  -  ,,  . .  —  A . ■  3  — 


o„  -0 

e  z 


o  -a 

z  r 


*Tre  2t9  z 


-  _  w 


or  in  terms  of  the  deviatoric  stresses  and  strains 


Yre  _  y8  z  _ 


yG  .  U) 


(8a) 


where  w  is  a  proportional  constant  that  depends  on  the 
position  of  the  volume  element  In  the  body  and  the  stage  of 
the  Inelastic  deformation.  For  «*1,  Hooke's  law  is  obtained. 
Octahedral  shearing  stress  t.  and  octahedral  shearing  strain 
y6  are  defined  by  ; 

tg-4‘  <9> 


(10) 
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The  third  hypothesis  gives  the  relation 


Ee  =  ( 1 -2p) S 


01) 


The  Prandtl-Reuss '  stress-strain  relations  embody 
three  hypotheses : 1 0 

1.  Principal  axes  of  stress  and  plastic  strain  incre¬ 
ments  coincide. 

2.  Mohr's  circle  diagrams  of  stress  and  plastic  strain 
increments  are  similar  at  any  stage  in  the  inelastic 
deformations. 

3.  Volume  changes  are  elastic. 

The  first  two  hypotheses  are  satisfied  by  the  relation 


lf.z  Sz  llri 

Sz  2V  2lrz 


(12) 


in  which  double  primes  denote  inelastic  components  of  strains 
and  dx  is  a  positive  scalar  factor  of  proportionality  which 
depends  on  the  position  of  the  volume  element  in  the  body 
and  the  stage  of  the  inelastic  deformation. 

The  third  hypothesis  can  be  written  as 


Ede  =  ( 1 -2y )ds 


(13) 


The  basic  differences  between  the  two  theories  cited 
above  have  been  described  in  References  10  and  15.  Hencky’s 
and  Prandtl-Reuss'  stress-strain  relations  yield  the  same 
results  if  the  principal  axes  of  stress  remain  fixed  in  di¬ 
rection,  otherwise  the  results  of  the  incremental  theory  will 
be  more  reliable. 

Experiments16  have  shown  that  the  inelastic  components 
of  strain  do  not  contribute  to  volume  change.  Stated 
another  way,  the  volume  change  is  always  elastic.  Hence, 


Ede  *  (l-2y)ds 
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If  the  stresses  satisfy  the  above-mentioned  equation,  the 
material  is  said  to  be  compressible.  A  material  is  said 
to  be  incompressible  if 

de  =  0 

To  reduce  the  complexity  of  the  mathematical  work,  the 
elastic  volume  changes  are  neglected.  In  this  investigation, 
however,  the  compressibility  of  a  material  will  be  taken  into 
account. 

THEORY 

1.  Assumptions  and  General  Information 

The  development  of  an  incremental  theory  for  thick-walled 
cylinders  is  based  on  the  following  assumptions: 

a.  The  material  is  homogeneous  and  isotropic. 

b.  The  deformation  is  small,  the  boundary  conditions 
at  the  surfaces  during  the  flow  are  determined  by 
the  radii  of  the  unstrained  cylinder. 

c.  The  inelastic  deformation  is  time-independent. 
d„  Prandtl -Reuss '  stress-strain  relation  is  valid. 

e.  The  material  is  assumed  to  be  compressible. 

f.  Von  Mises'  yield  criterion  is  valid. 

In  the  development  of  an  inelastic  theory  for  the  thick- 
walled  cylinder,  the  use  of  cylindrical  coordinates  (r,  e,  z) 
is  convenient  where  z  is  coincident  with  the  longitudinal  axis 
of  the  tube.  For  any  given  point  in  a  thick-walled  cylinder, 
the  nonvanishing  stress  components  are  of,  o0,  and  o2  while 

the  nonvanishing  strain  components  are  er,  eQ,  and  ez.  Because 

o *  radial  symmetry,  a  point  In  the  cylinder  has  only  two  com¬ 
ponents  of  displacement,  ur  and  u2,  except  for  rigid-body  dis¬ 
placements.  The  principal  strains  can  be  obtained  In  terms 
of  u^  and  u^,  and  are  given  by  the  relations: 


e 


r 


04) 


10 


(15) 


_r 

r 


3u . 

Jz 


(16) 


The  compatibility  equation  can  be  derived  from  Equations 
14  and  15  as  follows: 


de 

r”dr  =  Ve9 


(17) 


in  which  the  total  derivative  has  replaced  the  partial  de¬ 
rivative  since  er  and  eQ  are  independent  of  6  and  z. 

The  equation  of  equilibrium  may  be  written  as 


da, 

dr 


o„  -a 
0  r 


(18) 


The  differential  Equations  17  and  18  are  valid  for  both 
the  elastic  and  the  inelastic  regions  of  the  thick-walled 
cyl i nder. 

2.  Incremental  Theory  for  Elastoplastic  thick-walled  cylinders 

A  section  of  a  thick-walled  cylinder  is  shown  in  Figure  1 
in  which  the  plastic  region  has  spread  to  a  radius  c.  The 
cross  section  of  a  thick-walled  cylinder  is  divided  into  N 
rings  by  radii  r^a,  r,,  •  ••  rN+1«b. 

The  loads  applied  to  the  thick-walled  cylinder  will  be 
increased  step  by  step.  For  each  increment  of  load,  the 
finite-difference  method  will  be  used  to  solve  the  following 
formulated  problem: 

The  stresses  and  the  strains  at  any  point  r*ri  in  the 
elastic  region  must  satisfy  Hooke's  laws: 

nr  ■  iH°0  *  ^°z  *  d€r  *  er  -  r  °i))  (19) 


ii 


daA 

4  _ i 

Euvr  E 


U  0  ^ 

r  -  & a2  "  dee  ee  "  F  cVM(az  +  ar)] 


(20) 


fdar  ■  rda8  + 


■Z  "  r  LOe 


[o0-u(o^  +  ac ) ]  +  de_  (21 ) 


Prandtl -Reuss '  incremental  stress-strain  relations  are 
assumed  to  be  valid  at  any  point  in  the  plastic  region.  If 
the  volure  changes  are  assumed  to  be  elastic,  then  Equation 
12  reduces  to 


dej:  deg  de^ 


(22) 


Since  the  total  strain  in  an  overstrained  material  is  composed 
of  an  elastic  and  of  a  plastic  component,  i.e., 


e  =  e'  +  e" 


where  e'  and  e"  are  elastic  and  plastic  strain-components, 
respecti vely.  Then  Equation  22  can  be  rewritten  as 

<1(^-0  <*(e9-e^)  d(c,-e;) 


r  r 
S 


z  z 


(23) 


The  elastic-strain  components  e',  e ' »  and  e'  can  be  elimin- 

r  o  z 

ated  by  Hooke's  laws,  the  definitions  of  Sr*cr-S  * 

|-(2or-og  -c2) ,  etc.;  then  Equation  23  reduces  to 


der-  fkor-u(ce*o2)]  dce-|-[oe-u(oz+or)]  u(of+oe  )] 


?or-oe-o2 


?o0-  V°r 


W°e 
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Two  independent  equations  are  derived  as  follows: 


T  [<2VVV+“(2VVffz,ldVtt(2VVV'(2VV0z)1  d% 

£■  t(2dr-'3e-oz)+u(2o2-or-<Je  )]do2  +  (2oz-or-o9  )]dcr 

(2or"ae"aZ)  dc2  (24) 

and 


^C(2az-or-oe)-(2oe-oz-ar)]dor4[(2o2-or-<Je)+M(2o6-o2-or.)]dcie 

j*M(2VV°e)  +  (2VV0r)]dV(2VV°e)d£e 

(2ae ■°2_0r)dez 


The  loading  function  specifies  the  condition  under  which 
plastic  deformation  will  take  place  in  the  member.  Hill'-5 
proposed  a  loading  function  for  incremental  theories  which 
is  valid  for  all  states  of  stress.  For  linear  strai n-harden- 
ing  material,  the  loading  function  can  be  written 


a  *  ( i -a )o „  +  o£  e 

e 


(26) 


where  at  is  the  slope  of  the  straight  line  in  the  plastic 
region  and  a  may  be  considered  as  a  strain-hardening  factor 
for  the  material.  Since  t  ■  Z'  ♦  e" ,  Equation  26  reduces 
to  the  relation. 10 


o 


(27) 


in  which  a  is  the  effective  stress  and  is  defined  by 


0  *  ^  /(vVI  +  (V°;!)l  +  (V0r),+6(Trel+T9z,+Tzr,)  <28> 


and  e"  is  the  effective  inelastic  strain  and  is  defined  by 


e“  *  4  /(^-e5)2  +  (eJ.£p2+(e”-Ep2+|(Y"e2tY;22+YjrM  (29) 


The  inelastic-strain  components  can  be  eliminated  by  the 
relation 


e[!.  =  er“er  =  er*  f  [or‘U(oe+az)]  (30) 


and  two  similar  equations  for  e'l  and  e" 


Substituting  Equations  28,  29,  and  30  into  Equation  27, 
and  taking  the  derivative,  one  obtains  the  following 
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-  uw/ 

{~;(2ar-o0-oz)  +  (l+yH[(2er-e0-e2)-(i+ii)(2or-a0-cz)]}  ~T~ 


i  a 

+  {  “ (2aQ-az-ar)  +  (1+y)«J»[{2ee-er-e2)-(l+p)(2aQ-oz-or)]}  -£■“ 

2a 

i  da. 

+  {-r(2a2-ar-aQ)  +  (l+y)(J>C(2ez-fcr-ee)-(l+n)(2a2-or-ae)]>  -y- 
2a 


-  <J>C(2er-eQ-ez)-(1+y)(2ar-ae-az)]  der 


-  <j>[(2ee-ez-er)-(l+y){2a0-az-ar)]  deQ 


=  <|>[(2ez-Er-e0)-(l+}i)(2az-ar-a6)]  dez  (31) 


in  which 

=  i  _2_  1  (32) 

9  1 -a  e" 


The  equation  of  compatibility  (17)  and  the  equation  of 
equilibrium  (18)  are  valid  for  both  the  elastic  and  the  in¬ 
elastic  regions  of  a  thick-walled  cylinder.  The  finite- 
difference  forms  of  these  two  equations  are  given  by 


<ri  +  l-2lr'l  )<dar)i-<r1  +  i-r1 )  (dae>1  +  ri!d0'r)1  +  > 


*  <V0r)1(r1-M-r1)-r1[(°r)1+r<0r>lJ  (33) 
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for  the  equation  or  equilibrium,  and 


•(r'i+rri)(der)i+(ri+rri)(de0)i+ri(de0)i+i 


^£r"e6^i^ri  +  1"ri^”r’i^£e^i  +  i”^e0^i^ 


(34) 


for  the  equation  of  compatibility. 

At  each  point  r  =  r . ,  six  incremental  quantities  dar, 

danJ  da,,  de  .  deQ,  and  de,  exist  that  have  to  be  deter- 
u  z  r  y  z 

mined  for  each  step  of  increment  of  loading.  With  use  of  the 
fact  that  the  axial  strain  e  is  independent  of  r  if  the 
increment  of  axial  strain,  z  de  ,  is  specified  in  each  in¬ 
crement  of  load,  then  only  five  z  incremental  unknowns  exist 
at  each  point.  Hence,  5(N+1)  unknowns  exist  that  have  to  be 
determined  for  each  increment  of  loading.  The  five  equations 
listed  above  can  be  set  upat  each  station  (except  at  r=b), 
either  in  the  elastic  region  or  in  the  plastic  region.  At  the 
outer  surface  of  a  thick-walled  cylinder,  r=b,  Equations  33 
and  34  do  not  apply.  Hence,  the  total  number  of  equations 
is  5(N+l)-2.  For  the  solving  of  5 ( N+l )  unknowns,  two  ad¬ 
ditional  equations  resulting  from  the  boundary  conditions  are 


<di,r>r=ai 


<d0A.b! 


The  method  of  step-by-step  calculation  proceeds  as 
follows:  The  analysis  starts  with  known  loads  and  the  load¬ 
ing  path  is  divided  into  a  number  of  increments.  At  the 
beginning  of  each  increment  of  loading,  the  distribution  of 
stresses  and  strains  is  assumed  to  be  known.  To  save  com¬ 
puting  time,  the  stress  and  strain  distribution  for  a  thick- 
walled  cylinder  may  be  obtained  from  the  elastic  solution 
until  initial  yielding  occurs  at  the  inner  surface. 
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AP.j  and 


Step  1.  Specify  the  values  of  (doL  =  - 

r  r -a 

<^°r>r=b”  '  iPo' 

Step  2.  Assume  a  value  for  dez  (independent  of  r). 

Step  3.  Calculate  (dor)^,  {daQ)^,  (daz)..,  (def).j,  and 

(deQ)i  (for  1  =  1 ,2j*  •  *,N+1 )  from  a  5 ( N+l  )X5(N+1 ) 

matrix  which  is  formed  by  using  equations  (19),  (20), 
(21),  (24),  (25),  (31),  (33),  (34),  (35),  and  (36). 

Step  4.  Calculate** 

(tfyJi  -  (ar).j  before  an 
i ncreasing 
of  load 

(crQ ) i ,  (az)it  ( e r ) ^- ,  (eQ)i  and  (ez).  are  calcu¬ 
lated  in  the  same  way. 

Step  5.  Calculate  axial  load,  p  ,  by  using  Simpson's  law 

and  compare  with  actual  applied  axial  load  P  .  If 

a 

the  difference  of  P„  and  P,  falls  within  certain 

c  a 

allowable  limits,  the  calculated  stresses  and  strains 
in  Step  4  are  considered  to  be  acceptable;  then  go 
back  to  Step  1.  Another  set  of  new  increments  of 
pressure  and  strain  will  be  assigned  to  calculate 
a  new  stress  and  strain  distribution. 

If  the  difference  of  P_  and  Pa  is  greater  than  some 
allowable  limit,  a  new  dez  has  to  be  assumed;  then 
go  back  to  Step  2. 

The  elastic-plastic  boundary  is  determined  by 
o=  /(or-oe)2+(ae-az)2+(az-ar)2  »  og  (37) 

vc 


Step  6. 


Step  7. 
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In  the  elastic-plastic  problem,  the  elastic-plastic 
boundary  was  assumed  to  be  moved  In  an  increment 
equal  to  or  less  than  the  thickness  of  a  volume 
element. 


It  will  be  noted  that  a  system  of  equations  with  the 
nonzero  terms  clustered  about  the  main  diagonal  was  obtained. 
This  type  of  matrix  is  known  as  a  band  matrix  and  can  be 
solved  quite  rapidly  on  a  digital  computer.  In  the  computer 
program  which  was  developed,  the  Gaussian  elimination  method 
was  used  to  solve  these  equations.  All  calculations  were 
carried  out  with  double  precision  to  reduce  round-off  errors. 

EXAMPLES 

1.  Plane  Strain  Problem.  Most  of  the  theories  for  the 
elastoplastic,  thick-walled  cylinder  mentioned  in  the  intro¬ 
duction  are  derived  based  on  the  assumption  of  plane  strain, 
ez=0.  For  comparison  of  the  results  obtained  by  use  of  the 

present  theory  with  the  solutions  obtained  by  other  authors, 
the  problem  of  plane  strain  will  be  solved  first  in  this  in¬ 
vestigation.  For  the  purpose  of  obtaining  the  solution  to 
a  particular  problem,  the  following  numerical  values  are 
assigned  to  the  parameters  involved:  b/aa2,  p=0.3,  a=0.05, 
and  PQ=0.  The  cross  section  of  the  thick-walled  cylinder  is 

divided  into  20  rings  by  the  dimensionless  radii  p^=l, 

o2*1.05t"'p21«2.  A  105  by  105  matrix  for  each  Increment  of 

internal  pressure  (dar )p?|“Apl  *s  set  according  to  Equations 

19  through  36.  The  solution  to  this  particular  problem  can 
readily  be  obtained  by  use  of  the  procedure  stated  in  the 
previous  section  in  which  de2  Is  set  to  be  zero  for  each  in¬ 
crement  of  Internal  pressure.  The  axial,  radial,  and  cir¬ 
cumferential  stress  distributions  for  a  number  of  depths  of 
yielding  are  shown  In  Figure  3.  The  radial  and  circumferen¬ 
tial  strains  of  bore  surface  vs.  Internal  pressure  are  shown 
in  Figure  4. 

2-  Closed-End,  Th4ck-Walled  Cylinders.  A  closed-end , thick- 
walled  cylinder  is  a  cylinder  for  which  the  external  axial 
load  (exclude  the  axial  load  due  to  Internal  and  external 
pressure)  Is  zero.  To  obtain  a  numerical  solution  to  a 
closed-end,  thick-walled  cylinder  problem,  the  following 
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FIGURE  3.  Distribution  of  oxiol, radial,  and  circumferential 
stress  components  .of  .a  partially  plastic ,  thick-walled 
cylinder  subjected  to  internal  pressure 
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FIGURE  4.  Bore  Radial  and  Circumferential  Strains  of  a  Partly  Plastic  Thick-Walled 
Cylinder  Subjected  to  Internal  Pressure 


numerical  values  are  assigned:  b/a=l,5,  y=0,3,  a~0.05, 

and  P  =0.  The  solution  can  be  obtained  by  use  of  the  method 
o 

in  the  previous  section.  It  will  be  noted  that  for  a  closed- 
end,  thick-walled  cylinder,  the  only  axial  load  acting  on 
the  cylinder  is  due  to  the  internal  pressure  acting  on  the 
ends  of  the  cylinder.  Hence  dsz  will  be  found  such  that  the 

actual  calculated  axial  load  must  be  equal  to  (or  within  a 
certain  allowable  limit)  the  axial  load  due  to  internal  pres¬ 
sure,  The  distribution  of  axial,  radial,  and  circumferential 
stresses  is  shown  in  Figure  5.  The  axial  strain,  radial 
strain  of  bore,  and  circumferential  strains  of  bore  and  outer 
surfaces  of  a  thick-walled  cylinder  are  shown  in  Figure  6, 

3.  Open-End,  Thick-Walled  Cylinders.  A  solution  for  the 
open-end  cylinder  can  be  obtained  in  the  same  way  as  that 
for  the  closed-end  cylinder,  the  only  difference  being  that, 
for  each  increment  of  internal  pressure,  AP^  ,  the  calculated 

axial  load,  P  ,  has  to  be  zero  or  within  certain  limits  of 

V 

zero,  The  distribution  of  radial  and  circumferential  stresses 
is  shown  in  Figure  7.  The  axial,  radial,  and  circumferential 
strains  of  bore  surface  are  shown  in  Figure  8. 

4.  Closed-End,  Thick-Walled  Cylinder  Subjected  to  an  Addi¬ 
tional  Axial  Load.  The  purpose  of  this  example  is  intended 
to  show  that  the  solution  to  the  problem  of  an  el astopl asti c , 
thick-walled  cylinder  subjected  to  any  combination  of  in¬ 
ternal  pressure  and  axial  load  can  also  be  obtained  by  the 
present  method  without  any  change  of  digital  computer  program. 
The  problem  considered  here  is  that  of  internal  pressure 
applied  to  a  closed-end,  thick-walled  cylinder  until  yield 

Is  impending  at  the  bore  surface,  and  then  the  internal  pres¬ 
sure  Is  kept  constant  while  additional  axial  loads  are  ap¬ 
plied  at  both  ends  of  the  cylinder  to  move  the  elastic- 
plastic  boundary  from  the  Inner  surface  to  the  outer  surface. 
The  solution  to  this  problem  can  be  readily  obtained  by  use 
of  the  procedure  stated  In  the  previous  section.  The  dis¬ 
tribution  of  radial,  axial,  and  circumferential  stresses 
at  a  number  of  elastic-plastic  boundaries  Is  shown  in  Figures 
9  and  10.  Axial  load  F  vs.  axial  strain  and  circumferential 
strain,  and  radial  and  circumferential  strain  at  bore  surface 
are  shown  In  Figure  11.  (Axial  load  F  does  not  include  the 
axial  load  due  to  Internal  pressure.) 
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FIGURE  5.  Distribution  of  Axial,  Radial,  and  Circumferentl 
Stresses  of  a  Partly  Plastic,  Closed-End,  Thick- 
Walled  Cylinder  (y«0.3,  a»0.05,  a/b*T.5). 

2* 


1  a  1 ,  and  Circumferential  Strains  of  a  Partly  Plastic,  Ciosed-End, 
ed  Cylinder  (v-0.3,  a«0.05,  b/a~i,5) 


FIGURE  8  Axial,  Radial,  and  Circumferential  Strains  of  a  Partly  Plastic  Open-Ended 
Thick-Walled  Cylinder  fu-0.3,  a-G  05,  b/a-1  5), 


FIGURE  9  Distribution  of  Axial  Stress  of  a  Partly  Plastic, 
Thick-Walled  Cylinder  Subjected  to  Internal  Pressure  Until 
Yielding  Occurred  at  Bore,  and  then  Additional  End  Loads 
Were  Added  (y=0.3,  a=0.1,  b/a=1.5). 
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FI6URE  10.  Distribution  of  Radial  and  Circumferential  Stress 
?!  *  P^stic,  Thick-Walled  Cylinder  Subjected  to 

Pr®ssure  Until  Yielding  Occurred  at  Bore  and  then 
Additional  End  Loads  Were  Added  (u=0.3,  a-o.l,  b/a>1.5). 


FIGURE  11.  Axial,  Radial,  and  Circumferential  Strains  of  a  Partly  Plastic,  Thick- 
Walled  Cylinder  Subjected  to  Internal  Pressure  Until  Yielding  Occurred  at  Bore, 
and  then  Additional  End  Loads  Were  Added  (u=0.3,  a=0.1,  b/a-1.5) 


COMPARISON  OF  RESULTS 


The  problem  of  the  inelastic  behavior  of  a  thick-walled 
cylinder  subjected  to  internal  pressure  and  axial  load  has 
been  studied  from  various  points  of  view.  The  main  differ¬ 
ences  in  the  study  of  this  problem  are  those  of  stress-strain 
relation  in  the  inelastic  region,  yield  criterion,  allowance 
of  compatibility,  strain-hardening  of  material,  and  the 
assumption  concerning  the  axial  boundary  conditions.  A  survey 
of  the  various  theories  is  given  in  the  introduction.  Because 
of  the  difficulty  involved  in  the  estimation  of  error  for 
each  theory,  a  simple  comparison  will  be  used  to  show  the 
magnitude  of  the  corrections  introduced  by  the  present  theory. 
On  the  basis  of  the  assumption  of  plane  strain,  a  comparison 
is  given  in  Figure  12  of  the  stresses  predicted  by  the  pres¬ 
ent  theory  with  those  predicted  by  the  theories  of  (A)  Hill, 
Lee,  and  Tupper,1*  (B)  Hodge  and  White5  (Prandtl -Reuss '  stress- 
strain  laws),  (C)  Hedge  and  White6  (Hencky's  stress-strain 
laws),  and  (D)  Nadai. 

Curve  A  was  obtained  by  Hill  ,  Lee,  and  Tupper,  and  was 
based  on  Prandtl -Reuss '  stress-strain  law  with  the  Tresca's 
yield  condition  and  with  zero  strai n-hardening  ( a=0 ) . 

The  results  of  Hodge  and  White  are  presented  in  Curve  B. 
These  results  were  obtained  by  use  of  Prandtl -Reuss 1  stress- 
strain  law  with  Von  Mises'  yield  condition,  without  consider¬ 
ation  of  the  strain-hardening  property  of  the  material.  The 
results  obtained  by  Hodge  and  White  are  shown  in  Curve  C, 
and  are  based  on  Hencky's  stress-strain  law  with  Von  Mises' 
yield  condition,  and  without  consideration  of  strain-harden¬ 
ing  property  of  a  material. 

The  results  of  Nadai  are  given  In  Curve  D.  These  re¬ 
sults  were  obtained  without  consideration  of  compressibility 
of  the  material  in  the  plastic  region  of  a  cylinder.  Nadai 's 
solution  leads  to  the  physically  unrealistic  discontinuity 
in  the  axial  stress  at  the  elastic-plastic  boundary. 

The  results  obtained  by  use  of  the  present  theory  are 
shown  in  Curve  E. 

From  the  comparison,  good  agreement  between  predicting 
axial,  radial,  and  circumferential  stress  by  use  of  the 
present  theory  and  the  theories  of  (A),  (B)  and  (C)  has  been 
found . 
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FIGURE  12.  Comparison  of  Distribution  o*  Axial,  Radial,  and 
Circumferential  Stresses  Predicted  by  the  Theories  of  (A) 

Hill,  Lee,  and  Tupper,  (B)  Hodge  and  White  (Prandtl -Reuss 
Stress-Strain  Law),  (C)  Hodge  and  White  (Hencky  Stress-Strain 
Law),  (D)  Nadal,  and  (E)  Present  Approach  for  a  Partly  Plastic 
Thick-Walled  Cylinder  Under  Plane-Strain  Condition. 
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The  stresses  for  a  closed-end  cylinder,  predicted  by 
the  present  theory,  are  compared  in  Figure  13  with  those 
predicted  by  Smith  and  Sidebottom' s 1 0  theory.  The  greatest 
difference  is  noticed  in  the  axial  stress  at  the  bore  sur¬ 
face.  The  radial  and  circumferential  stresses  agree  fairly 
well  with  these  two  theories.  The  axial  and  the  circum¬ 
ferential  strains  at  the  outer  surface  of  a  cylinder  pre¬ 
dicted  by  the  present  theory, and  Smith  and  Sidebottom's 
theory  are  compared  in  Figure  6.  The  difference  in  the 
stresses  and  the  strains  in  this  comparison  is  obtained  by 
use  of  Prandtl -Reuss '  and  Hencky's  stress-strain  relations. 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  closed-form  elastoplastic  solutions  of  thick-walled 
cylinders  can  only  be  obtained  based  on  Hencky's  total  stress- 
strain  relations.  On  the  basis  of  Hencky's  laws,  the  strains 
and  stresses  at  any  state  are  dependent  upon  the  Instantaneous 
state  of  stresses  and  strains  and  are  not  dependent  upon  how 
that  stress  and  strain  system  is  reached.  Taylor  pointed  out 
the  inaccuracy  of  the  total-strain  assumption,  particularly 
for  a  member  subjected  to  nonproportionate  loading,  A  more 
rigorous  general  solution  to  the  problem  of  an  elastoplastic 
thick-walled  cylinder  subjected  to  either  proportionate  or 
nonproportionate  loading  can  readily  be  obtained  numerically 
with  comparitively  few  restrictive  assumptions  by  use  of 
the  theory  developed  in  this  report. 

The  finite-difference  treatment  of  the  elastoplastic 
problem  of  a  thick-walled  cylinder  is  normally  based  on  the 
differential  equations  for  the  displacement  vector;  there¬ 
fore,  determination  of  stresses  and  strains  requires  numerical 
differentiation.  However,  good  results  in  differentiation 
are  not  provided  by  the  computer  unless  a  rather  fine  grid  is 
used.  With  the  alternative  method,  developed  in  this  report, 
incremental  stresses  and  strains  at  each  nodal  point  are 
directly  used  as  variables;  hence,  numerical  differentiation 
in  the  evaluation  of  stresses  and  strains  is  not  required. 

The  theory,  presented  In  this  report,  includes  Prandtl- 
Reuss'  incremental  stress-strain  laws,  Von  Mlses'  flow  cri¬ 
terion,  and  the  strain-hardening  and  compressibility  properties 
of  the  material.  The  stresses  and  the  strains  In  all  principal 
directions  can  be  computed  at  the  same  time;  discontinuities 
of  stresses  or  strains  are  not  presented  at  the  elastic- 
plastic  boundary. 

The  comparison  of  distribution  of  stresses  predicted  by 
the  present  theory  with  those  predicted  by  the  theories  of 
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Hill,  Hodge,  and  Nadai  was  made  under  the  plain-strain  con¬ 
dition  without  consideration  of  strain-hardening.  Excellent 
agreement  was  found  amoung  the  stresses  predicted  by  the 
theories  of  Hill,  Hodge,  and  the  present  theory.  However, 
Nadai's  solution  leads  to  physically  unrealistic  discontinuity 
in  the  axial  stress  at  the  elastic-plastic  boundary. 

For  the  expansion  of  the  application  of  the  developed 
theory  to  actual  engineering  design  problems,  the  following 
recommendations  are  suggested: 

1.  Experimental  evaluation  of  the  theory  should  be 
performed . 

2.  The  inelastic  behavior  of  thick-walled  cylinders 
subjected  to  cyclic  loading  should  be  Investigated. 

3.  The  present  theory  should  be  extended  to  include  the 
effect  of  thermal  loading. 


FIGURE  13  Comparison  of  Distribution  of  Axial,  Radial,  and 
Circumferential  Stresses  Predicted  by  the  Theories  of  (A) 
Smith  and  Sidebottom  (Hencky  Stress-Strain  Law);  (B)  Present 
Approach  (Prandtl -Reuss  Stress-Strain  Law)  for  a  Partly 
Plastic,  Closed-End,  Thick-Walled  Cylinder  (w-0  3,  a*0.l, 
b/a  =  1 . 5) . 
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